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Abstract 

The paper aims to generalize Clausen's identity to the square of any Gauss 
hypergeometric function. Accordingly, solutions of the related 3rd order linear 
differential equation are found in terms of certain bivariate series that can reduce 
to 3F2 series similar to those in Clausen's identity. The general contiguous variation 
of Clausen's identity is found. The related Chaundy's identity is generalized without 
any restriction on the parameters of Gauss hypergeometric function. The special 
case of dihedral Gauss hypergeometric functions is underscored. 
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1 Introduction 

Clausen's identity [Cla28| is 



2F1 



a, b 
a + b + 



2a, 2b, a + b 
2a + 2b,a + b + 



A related identity is |Cha58[ (5)]: 



2F1 



a, b 

a + b + 



2F1I 



2 

a - 




(1) 



(2) 



Up to a power factor, the two 2F1 functions in the second identity are solutions of the 
same Euler's hypergeometric differential equation. The two identities demonstrate a case 
when the symmetric tensor square of Euler's hypergeometric equation coincides with the 
hypergeometric differential equation for 3F2 functions. Accordingly, formulas ^ and ^ 
relate quadratic forms in 2F1 functions and 3F2 functions as solutions of the same third 
order linear differential equation; see Exercise 13 in [AAR99[ pg. 116]. 

This paper aims to generalize Clausen's identity to the square of a general 2F1 func- 
tion, without any restriction on the three parameters. Accordingly, we try to find other 
attractive solutions of the symmetric tensor square equation for general Euler's hyper- 
geometric equation. We present the symmetric tensor square equation in Section [3l see 
formula (ISTI). 
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We find that the general symmetric square equation has solutions expressible as spe- 
cializations of the following double hypergeometric series: 



2:1:1 / a;b;pi,p2 



1:1;1 



c; 91,92 



Fr 



2:0:0 



a;pi,P2;gi,g2 

b\ c 



x,y 



x,y 



oo oo 



t'oh (c).+.(9i).fe),z!j! ^' 



oc> oo 

EE 

i=0 3=0 



{a)i+j jpi), {p2)j {ql)^ {q2)j 



(3) 
(4) 



li a = c, these two series reduce to the bivariate Appell's F2 or F3 hypergeometric series, 
respectively. 

Surely, the left-hand sides of (IT])-(l2]) are trivially double hypergeometric series as 
well. The significance of new solutions is that they reduce to 3F2 functions of similar 
shape as in ([T])-([2]) in special cases, and that they represent elementary solutions (that 
is, a polynomial times a power function) if the the symmetric square equation has such 
solutions. 

Our main results are summarized in the following two theorems. In Theorem II. 1[ 



by a univariate specialization of F^.'^.'^ or F^'.^!^ function we understand the restriction 



2:1:1 



of any branch (under analytic continuation) of the bivariate F^.'^J^ or F^'.^'^ function to 
the curve parameterized by z. In Theorem II. 2i we refer to the F^ll^l or F^.'^J^ series as 
defined in ([3|)-(l4|). A direct generalization of Clausen's identity is formula (fTO| . It applies 
to the Gauss hypergeometric functions contiguous to Clausen's instance. We prove these 
theorems in Sections S] and [5l respectively. 

Theorem 1.1 The univariate functions 



F, 



1:1;1 



= (l-z)' 



;_a-b-ipl:2;2 
-^2:0:0 



2a; 2b; c - ^, a + b- c + ^ 
^ a + 5 + i; 2c - 1, 2a + 2& - 2c -I- 1 

-6+i,a+6— c+i;6— a+i,c— a — 



= (l-z) 



':—a — b T7l:2;2 



^2:0:0 



c; 2 - C 

^; a— 6+-i,c— i; 6 — a 
a + fe — c+1; c — a — 



z-l 



' 2 ' 2 
6+1 



l-z,l-- 
z , 



z^ 



-1-6 — c+i;| — c, c— a 
a — 6+1: 6— a + 1 



6+i 



1 1 



z 1 — z 



(5) 
(6) 
(7) 
(8) 



satisfy the symmetric tensor square equation for 2Pi{^'c ■ 
Theorem 1.2 The following identities hold in a neighborhood of z — 0; 
1 + a — c, 1 + 6— c 



a, 6 
c 



2:0:0 



2-c 
6 + 5, a 



= {I - z)"-"-^—^ X 



c+i;6 
c; 2 - c 



6+i 



z-l 



(9) 



2 



2Fl 



a, b 

a + b + n + ^ 



{^)n [o- + b + j-,2:l;l 



(a+ i)„ (6+ i), 



■-^1:1:1 



(a + ' 



2-ri 



a, —a — m — 71 
1 ^ 



2a; 26; a + 6 + n, 
2 



i; 2a + 2& + 2n,-2n 



(^) 



Fi 



1 1 

2' 2 



m+1 V 2 / n v2/m+n 



1:1;1 



2a; —2a — 2m — 2n; 



-m, 



— 771 — 71 ; — 277T,, —2n 



z, \ — z 



(10) 



(11) 



Here a, 6, c can be any complex numbers, if only the involved lower parameters c, 2 — c, 
a + 6 + ^ are not zero or negative integers. The numbers m, n are assumed to be non- 
negative integers. 

The F^: i.\ series in (fTO|) is understood to terminate in the second argument at the power 
(1 — z)", but it is never a terminating series in the first argument even if a + 6 + 7i is a 
non-positive integer —771 (when the term- wise limit b —a — n — m should be applied). 
The i^i^.i.'i series in (|lip is a finite rectangular sum with [m + l)(n + 1) terms. 

If 71 = 0, formula (fTO|) reduces to Clausen's identity ([T]). Ifc = a + fe+^' formula 
© reduces to ^ . We should observe that Clausen's identity is wrong if a -I- 6 is a non- 
positive integer and the 3F2 series is interpreted as terminating. A correct identity with 
the terminating 3F2 series is then the special case 77 = of pT|) : 



2F1 



3^2 



2a, 



\2)m \ 2)m+l 

-2a — 2m, —m 

— 771, — 2?7l 



2F1 



777, -f 
777, 



1 1 

2 ' 2 

^ 2 



(12) 



Formula (fTTj) has the following significance. The 2F1 functions on the left-hand side 
are dihedral Gauss hypergeometric functions, as the monodromy group of their Euler's 
hypergeometric equation is a dihedral group. Then the symmetric tensor square equa- 
tion is reducible, and its monodromy representation has an invariant one-dimensional 
subspace. The F'^:^'^ function on the right-hand side is a polynomial, so it is obviously 
an invariant of the monodromy group. Formula (|lip identifies a generator for the in- 
variant space as a linear combination of 2Fi(z)^ solutions, and as an explicit terminating 
Fi.'i.'i expression |Vid08|, Section 5]. In more plain terms, the Fj^.i.'i expression gives 
an expected elementary solution of the symmetric square equation. Similarly, the F2.q.'q 
series in ([9]) is terminating in both summation directions in the following specialization 
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with dihedral 2F1 functions: 



2-ri 



a, a + 771 + i 

2a + 777 



2 

77+1 



2F1 



-a — 777 — 77 

1 - 2a 



' 2 
777 



a - 

77 



Z ] ^{l- Z)" X 



^1:2;2 



2a 



— 777, —77; 777 + 1, 77 + 1 
77 + 77 +1:1 — 2a — 777 - 



z, ■ 



(13) 



The terminating F^.'^.'^ and F2.Q.Q sums are closely related. Up to a constant multiple, 
one can rewrite the terminating F'^'^'^ sum as a terminating F2.'q.'q sum, and vice versa, 
simply by reversing the order of summation in both directions. 



2 Variation of formulas 

A product of two single hypergeometric series is trivially a double hypergeometric series. 
In particular, we can write 



2F1 



a, b 



n0:2;2 / O.] 6, b 

c,c 



Chaundy |Cha58[ (6)] gives the following double series expansion 

2 oo 



2F1 



a, b 



^ (2a)fc(26)fc(c-l)fc /-|,-fe^,i,a+&-c+i 
(c)fc(2c-l)fcfc! ^ a 



1 1 -r- k 

2 ' 2 



1 z" 



(14) 



(15) 



But our formulas (dU]) and ^ are particularly interesting as direct generalizations of the 
classical formulas ([I])-©. 

Application of Euler's transformation [AAR99[ (2.2.7)] to the second 2F1 factor on 
the left-hand side of (O gives 



a, b 



z 2F1 



1 - a, 1 - 

2-c 



pl:2;2f 2 
2:0:0 



i;a — 6+i,a + 5— c+i;6 — a+i 



2, c — a — 6+2 



c; 2 - c 



z, 



z-l 



(16) 



The substitution z i-^ z/(z — 1) and application of Pfaff's transformation [AAR99[ (2.2.6)] 
to the 2F1 functions in pT|) and (fT^ gives, respectively: 



(a+i)Ja + 777 +i) 



' n 



2y„ / a, a + n+ 2 

21*1 1 

\ 7^—777 



(«+2)rn(° + "+l)m 2m+l ^ f a + 777 +i,a + 77 + 777+l 



(2)™ (2)m + l (2)11 (2)m+„ 



(1-z) 



-2a p2 



■ i-i 2a; —2a — 2777 — 277; —777, —77 

■i — 777 — 77; —2777, —277 



z- 1' 1 -z. 



(17) 
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2Fl| 



a, a + ^ 



(2)™ (2) 



- 2^1 



Notice that 



(1-Z)-^«3F2 



2a, —2a — 2m, — m 



i — TO, — 2to 



1 



T2m i 



(to!)2 (2a)2,„^ 



^2)7™ (2)r?x+l 



(2m)!(2m + l)! 

The speciahzations to = of and (fT7|) are interesting as well: 

2 



(18) 
(19) 



2yn 



2F1 



a, —a — n 
1/2 

'1 n2 - 



3J:^2 



a, a + n + i 
1/2 



2a, —2a — 2n, —n 
i — n, — 2n 
2 

2 



2z(a)^ 



>2 s_2a " / 2a, 



(^):(i--r"3F2 



r("-^^ 


,i-a- 




3/2 








- i, a + 


2Fr(«^ 


3/2 


2n, — n 




-2n 





(20) 



(21) 



It is tempting to equate 2Fi^''^^ to the F^.'^'l fmiction in ([5]) up to a constant 
multiple, since the latter appears to have the power series expansion at z = as well. 
The fcth derivative of the F^.'^'^ function at z = would evaluate to a linear combination 
of the values 



3^2 



1 



j=0,l,...,/c. 



(22) 



2a + k, 2b + k, a + b - c + j + ^ 
a + b + k+^,2a + 2b-2c + j + l 

However, the convergence condition of these 3F2(1) sums is Re(l — c — fc) > 0, so the 
3F2(1) values are undefined for large enough fc. The ff.'i.'i function can have branching 
behavior at z = in general, as the local exponents z^~'^, z^"^"^ of the symmetric square 
equation can come into play. See [Vid09a| for explicit details. 

The Fi/i.'i function in © can be evaluated at z = and z = 1 using the 3F2(1) 
evaluation in ' AAR99[ Theorem 



1:1:1 



2a; 26; c 
-b 



i,a + 6-c+i 



0,1 



p2:l;l| 
-^1:1:1 



i; 2c- l,2a + 26- 2c+ 1 , 

r (^) r (a + 6 + i) r(i - c) r(i + a + b - c) 
r (a + i) r (6 + i) r(i + a - c) r(i + 5 - c) ^ 

2a; 26; c - i, a + 6- c ' ^ 
- 6 + i; 2c - 1, 2a + 26 - 2c + 1 

r (I) r (a + 6 + ^) r(c) r(c - a - 6) 
r (a + i) r (6 + i) r(c - a) r(c - 6) ' 




(23) 



(24) 
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The convergence conditions are Re(l — c) > and Re(c — a — b) > 0, respectively. 
It is worth mentioning that Bailey's identity [Bai33] 



F4 

specializes to 



a; b 
c, a+b—c+1 



- 2/),y(l - x) 



a; b 
c, a+b—c+1 



2F1 



2F1 



a, b 
c 



a, b 
c 



a; 2F1 



a, b 
a+b-c+1 



2: 2F1 



a, b 
a+b~c+l 



1-x 



(25) 



(26) 



The Fi^x^, (1 — x)^) function is another solution of the symmetric tensor square equation. 
In particular, we have the following identity for x S [0, 1] if Rec < 1 and Re(c— a — 6) > 0, 
due to the connection formula in [AAR99| (2.3.13)]: 



2^1 



a, b 
c 



T{c)T{c- a-b) 



Fa 



T{c-a)T{c-b) "Kca+b-c+l 
(1 - tY''-'"'-"'Fi 



x\{l~xY\ + 



T{c)T{a + b-c) ^,2c-2a-2b 



r(«)r(5) 

3 The differential equations 

Let 02, Qx, 01/ denote the differential operators 



c — a; c — 6 
c, c— a — 6+1 



x^,{i~xy 



(27) 



dz ' 



d_ 

dx' 



Qy = V 



dy' 



(28) 



Euler's hypergeometric differential equation [AAR991 (2.3.5)] for a general Gauss hyper- 

a, b 



geometric function 2F1 



can be compactly written as follows: 



z (6, + a) (6, + b) y{z) - 9, (9, + c - 1) y{z) = 0. 



(29) 



This is a Fuchsian equation with three regular singular points z = 0, 1 and 00. The local 
exponents at the singular points are: 

0, 1 — c at z = 0; 0, c — a — b at z — 1; and a, 6 at z = 00. 

In general, a basis of solutions is 

1 + a — c, 1 + 6— c 



a, b 



2-c 



(30) 



Note that the left-hand sides of formulas ^ and © are products of these two solutions 
of the same Euler's equation, with the power factor z^""^ ignored. The classical formulas 
([T])-© apply when the difference of local exponents at z = 1 is equal to 1/2. 
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If 2/1 , 2/2 are two independent solutions of a second order homogeneous linear differen- 
tial equation like (|29|) . the functions y1, yiy2, satisfy a third order linear differential 
equation called the symmetric tensor square of the second order equation. The symmet- 
ric tensor square equation for Euler's equation can be written as follows [Cha58i 
(2)]: 

z {Q, + 2a) (9, + 2b) (9, + a + b) y{z) - 949, + c - 1)(9, + 2c - 2) y{z) 

_^ {2a + 2b-2c+\)z = 0. (31) 

z — 1 ^ ' 

Ifc = a-|-5-|-i, this is a differential equation for the 3F2 function in ([T]). 

Partial differential equations for the F-^:-^!-^ (x, y) function in ([3]) can be obtained by 
considering the first order recurrence relations between its coefficients. Writing the fI:1:1 
sum as X^i^o X^^o '-'J ' have the relations 

Cj+i,] _ [a + i + j){b + i + j){pi +i)x Cjj+i _ {a + i + j){b + i + j){p2 + j) y 
Ci,3 {c + i + j){qi+i){l + i) ' c^j (c + i + j) (92 + i) (1 + j) 

also Q+i j/ Cij+i = (pi + i){q2 + + j) xj {p2 + j){qi + i)(l + i) y. These relations 
translate to the following partial differential equations for the F^.'i.i (x,y) function: 

Pi = a;(9, + 9y + a)(9, + 9y + 6)(9, +pi) 

-9,(9, + 9j, + c-l)(9, + (7i-l), (32) 
P2 = y(9, + 9y + a)(9, + 9y + 6)(9y+p2) 

-9j;(9, + 9j; + c-l) (9^ + 92-1), (33) 

P3 = xey{ey + q2-i)ie^+pi)-ye,ie, + qi~i){ey + p2). (34) 

In general, two of these operators generate the ideal J in C{x,y){Qx,Qy) annihilating 
the F^.'l'l {x , y) series. In particular, we have the following obvious syzygies: 

-2/(9y+p2)Pi+a;(9, +pi)P2 + (9, + 9.y + c-2)P3 = 0, 
ey(9y + (72-l)Pi-9,(9, + (7i-l)P2 + (9,+ 9y+a-l)(9,+ 9y+6-l)P3 = 0. 

If the coefficient 92, + 9.y + c — 2 does not divide (9a; + Qy + a ~ l){Qx + Qy + b ~ 1), 
we can express P3 in terms of Pi, P2, etc. 

Note the following commutation relations in C{x, y){ldx; 'dy)'- 

exx = xex + x^x{ex + i), eyy = y{ey + i). (35) 

The variables x, Qx commute with y, Oy. Up to the multiplication order, these relations 
are the same as shift operator relations, such as Snn = (n + l)Sn in the 9re algebra 
C{n,Sn)- We will keep working with Euler-type differentiation operators in (|28p in our 
transformations of differential equations. 
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Grobner basis computations show that third order partial differential operators in the 
ideal J are hnearly generated by Pi,P2,^3- A residue basis over C{x,y) is formed by 
the 7 monomials 1, 9a;, 9^, G^jOj,, 0^, Bj^. Hence the rank of the differential system 
for F'^'^.'^ {x,y) functions is equal to 7. The leading coefficients in various Grobner bases 
suggests that the following lines are in the singular locus of the differential system: 



X = 0, a; = 1, y = 0, y = 1, 



y = l, y + 2x = l, 2y + x = l. 



It is easy to apply the transformations F{x,y) i— > x^y^ F{x,y) or x i— > 1/x, y i— > 1/y 
to the differential system generated by the operators Pi,P2. In effect, we only have 
to replace additionally ^ + ct, Qy i-^ Qy + P or 9a; i— > —92,, 9y ^©yi 
respectively. We may attempt to find transformations of the operators Pi , P2 to a pair 
of similar hypergeometric operators. This gives us the following solutions of the same 
system of partial differential equations, beside the F'^:^.'^ {x,y) function in ([3]): 



1-51 p2:l;l/ 1 
^1:1:1 



,1-92 p2:l;l 
/ ^1:1:1 



- gi ; 1 + 5 - gi ; 1 + pi - gi , p2 
1 + c ~ qi; 2 - qi, q2 
2:1.1 f l + a-q2\ l + b- q2;pi,l+P2 - q2 
1 + c - 92; 91, 2 - q2 

l-q^ l-q^ p2:l;lf 2 + a-qi-q2; 2 + 6-^1-^2; l+Pl-^l, 1+^2-92 

^ '^'''V 2 + c-qi-q2;2-qi,2-q2 

1:2:2 f ^ + Pi + P2 - C; 1 + pi ~ qi , 1 + P2 - q2;Pl,P2 



~P1 



2:0:0 



I+P1+P2- a; 1+P1+P2 



x,y 

x,y 

x,y 
1 1 

x' y 



(36) 
(37) 
(38) 
(39) 



As we see, the systems of partial differential equations for the F^.'^'^ and F2.Q.Q functions 
are easily transformable to each other. In the other direction, here is a F^.'i.i function 
satisfying the same partial differential equations as the F2.Q.Q {x,y) function in ([4]): 



-pi ,,-P2 p2:l;l I 1 
y ^1:1:1 



P1+P2- a; l+pi 
1 + Pi + P2 - c; 1 + pi - 



^P2-b; pi,P2 
qi,l + P2 ~ q2 



1 1 

x' y 



(40) 



From here, the other F^.'l.'l companion solutions for the ^2.0.0 (2;, y) function in ^ can 
be obtained using the symmetries pi <-> qi and P2 ^ q2- 



4 Proof of Theorem 11.11 



First we prove that the F^:^\ function in (O satisfies the symmetric square equation ([21 



for2Fi^°^^ . The ideal of partial differential operators annihilating 



p2:l;l 
-^1:1:1 



2a; 2&; c - i, a + 6- c + i 



; 2c -1,2a + 26- 2c +1 



x,y 



(41) 
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is generated by the polynomials Pi, P2, P3 in ([S^ - fM]) with the parameters specialized 
to linear functions in a, 6, c. Here are the specialized differential equations: 

= a; (e^ + + 2a) (e^ + e,y + 26) (6:, + c - i) 

-e,(e, + e^ + a + 6-i)(e, + 2c-2), (42) 

^ y{Q^ + Qy + 2a){e^ + ey + 2b){ey + a + b-c+^) 

-Qy {e^ + Oy + a + b- i) (Gy + 2a + 26 " 2c), (43) 

= xey{Qy + 2a + 2b-2c){Q^ + c-^) 

~yeAQx + 2c~2)iey + a + b-c+^). (44) 



The variables and derivatives are related as follows under the univariate specialization 
under consideration: 

x = y = i-z, ^ = ^-^, e. = e. + ^e,. (45) 

Following )Vid09b|, Definition 1.1], the partial differential form of an ordinary differential 
equation under a specialization like ()45ll is the expression where the univariate derivatives 
are replaced by respective linear combinations of partial derivatives. On the other hand, 
the specialized form of a partial differential equation under the same kind of specialization 
is the expression with coefficients to the partial derivatives specialized to univariate 
functions. In our setting, the two forms are "mixed" linear differential expressions in 
the partial derivatives Q^, Qy but with the coefficients univariate in z. Algebraically, 
we work in the C[z]{Qz) module generated by the partial derivatives of any order. The 
action of 83 on the partial derivatives is given by the identification in (|45p. 

To show that the univariate F'^'^.'^ function in (jlOp satisfies the symmetric tensor 
square equation (|3ip . we demonstrate that the partial differential form of the symmetric 
square equation coincides with the specialized form of a partial differential equation fol- 
lowing from the partial differential equations In the terminology of [VidOQbi 
Definition 1.1], we show that the symmetric square equation follows fully from the dif- 
ferential equations ((42|) - (|44)) under the speciafization defined in ([45]) . 

To give a partial differential form of the symmetric square equation l|3ip . we intro- 
duce a way to work with multiplicative expressions of (non-commutative) differential 
operators. Suppose that G, H are functions in z. Then 



(9, + G) (9, + H) = (9, + G) [e^ + —J Qy + H 

8. + -A_e, + A).(9. + ^9, + i.)-^9, + .f, (46) 

where we use the dot between two differential factors to signify commutative multiplica- 
tion. In other words, an expanded expression for (8z -I- G) {Qz + H) can be obtained by 
multiplying commutatively the two large factors in (j46p . collecting terms to the mono- 
mials in 93;, Oy, then interpreting those monomials as differential operators of suitable 
order, and writing the coefficients to those monomials as left-side factors. 
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If A, B, C are constants, the product {Qz + A) (8z + B) (8z + C) can be written as 
follows: 



[z-l) 



2 



2z 
z- 1 



e„ + B + c 



+ 1) 
" (z-l) 



0. 

3 



and finally, as 



1 

z 



(z-l)2 



z(z + l ) 
(z-l) 



•(36:^ + ^^Qy + A + S + C 



(47) 



Following this "commutative" expression, we can write the symmetric square equation 
(1311) as follows: 



z - 1 



z 



- e 



z-l 



e„ • e 



z-l 
z 



-e„ + 26- 6, 



z- 1 



e„ + c - 1 • 



z-l 



-e„ + 2c- 2 



3z 



(z-l) 
(2a + 26 - 2c+ 1) z 
z- 1 



z- 1 



-©. 



(a + fe)z + l-c^ z(z + l) 



z- 1 



T©. 



(a + 6- c+ 1) ( 6:^ + ^ j-Oy ) + 2a6 ) 0. (48) 



Using formal commutative computations, one can check that this equation coincides with 
the specialized (under x = z, y = 1 — z) form of 



[ Eq. dm) ] - [ Eq. ill) ] - [ Eq. gl ]. 



(49) 



We proved that the F'^:^.'^ function in ([5]) satisfies the symmetric square equation PT|) . 

The relation between the functions in ^ and (PO)) implies that the function in ^ 
satisfies the same symmetric square equation. The symmetric square equation has the 
following two solutions as well: 



Z 2ri 



a, 1 + a — c 
l + a-b 



(l-z)-^'^2Fi 



a, c ~ b 
l + a-b 



1 - z 



(50) 



The relation between the functions 2Fi^°^'' z^ and ([5]) translates into the claim that 
([6|) and ([7]) are solutions of the same symmetric square equation as well. This completes 



the proof of Theorem 11.1 
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Notice that the upper parameters in ([SI)-® contain the local exponent differences of 
2Fi(^'^'J' increased by 1/2, while the lower parameters are the same local exponent 

differences increased by 1. A set of twelve F'^:^'^ solutions can be obtained using the 
relation between the functions in ([4]) and (|40)). 



5 Proof of Theorem 11.21 

To prove formula ([5]), we compare the function in © with the following solution of the 
symmetric square equation: 



z'-' 2F1 



a, b 
c 



2F1 



1 + a-c, 1 + 5- 



(51) 



Both solutions have the same local exponent 1 — c at z = 0. For those values of c for 
which both 2F1 solutions are generally defined, the linear space of solutions with this 
local exponent is one-dimensional. After division by z^~^ both solutions ([6]) and (|5T|) 
evaluate to 1 at z = 0, therefore formula ^ follows. 

To prove formula pil)) . we observe that (for general a,b) the linear space of power 
series solutions at z = to the symmetric square equation is one-dimensional. The 
^i-i'i function in (jlOp and the square function in (I14p are both power series solutions 
with the local exponent 0. Hence they must differ by a constant multiple. To compute the 
constant multiple, we evaluate both sides oi z — 1. Note that the convergence condition 
Re (n -I- i) > for both sides is always satisfied. The left-hand side is evaluated by the 
Gauss formula [AAR99i Theorem 2.2.2]; the right-hand side is evaluated by ([M)) with 
c = a + b + n+^. The Pochhammer-type factor follows. 

To prove formula (fTT|) . we take the term- wise limit b —a — m — n in (fTO|) . We use 



lim 



(e - k)2k+i 



{-k)kk\ {-lf{k\f {-lfk\ 



£-0 (2e- 2fc)2fc+i (-2fc)2fe-2 
to get the limiting formula 



2 • (2fc)! 



(52) 



(2)71 (2 



m — n 

:i2:l;l 
1:1:1 



2F1 



a — m 
1 
2 



2a; —2a — 2to — 2n; — m, ~n 



-2m, — 2n 



zA 



(2a)2,n+i (-2a -2m- 2n)2m+i (-1)" m\ 
(^---")2™+i(2- + l)! 22^+1(1), 
2a + 2to + 1; 1 - 2a - 2n; m + 1, - 



Ff 



-\- m — n; 2m + 2, —2n 



^2m+l 



Z, 1 



(53) 
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Due to (|10p . the latter F^.'^J^ function can be written as 

(a + m + 1),^ (1 — a — rt) 



(i)Jm-n+|)„ 



2F1 



We replace (2m + 1)! ~ 2^™+^m! (i)m+i and rewrite identity ([55)1 as follows: 



(a + m + 



2F1 



a, —a ~ m — n 
1 ™ 



,2:i;i / 2a; —2a — 2m — 2n; — m, — n 



-^1:1:1 



zA — z 



^ ~ m ~ n; —2m, — 2n 

(-1)"+"+^ (2a)2n^+l (2a + 2n)2^+i 2,„+i 
24™+^(i)™(i)™+i (i-m-n)^„^^^ 



(a + m + l)„(a)„ 



(i)Jm-n+|), 



2F1 



a + m+ i,i — a — n 



After regrouping 



(54) 



(55) 



(i-m-n)^ (m-n+l) =(i-m-n), = (-1)™+" (i) (i) 

(2a)2m+i = 22™+i (a),„+i (a + \) ^ , (a),„4.i (a + m + 1)„ = (a)„+„+i, etc., 

we obtain (fTT|) . 

Alternatively, we may find the constant multiple in (jlOp by evaluating both sides at 
z = 0. The evaluation of the right-hand side can be done by Zeilberger's algorithm. 
Formula (jlip can be idependently proved by observing that the linear space of power 
series solutions at z = to the symmetric square equation is two-dimensional, and 
evaluating at the two points z = and z = 1. 
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